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ABSTRACT 

The response of cloud simulations to turbulence parameterizations is studied systematically using the GISS 
general circulation model (GCM) E2 employed in the Intergovernmental Panel on Climate Change’s (IPCC) 

Fifth Assessment Report (AR5). Without the turbulence parameterization, the relative humidity (RH) and the 
low cloud cover peak unrealistically close to the surface; with the dry convection or with only the local tur- 
bulence parameterization, these two quantities improve their vertical structures, but the vertical transport of 
water vapor is still weak in the planetary boundary layers (PBLs); with both local and nonlocal turbulence 
parameterizations, the RH and low cloud cover have better vertical structures in all latitudes due to more 
significant vertical transport of water vapor in the PBL. The study also compares the cloud and radiation cli- 
matologies obtained from an experiment using a newer version of turbulence parameterization being developed 
at GISS with those obtained from the AR5 version. This newer scheme differs from the AR5 version in com- 
puting nonlocal transports, turbulent length scale, and PBL height and shows significant improvements in cloud 
and radiation simulations, especially over the subtropical eastern oceans and the southern oceans. The diagnosed 
PBL heights appear to correlate well with the low cloud distribution over oceans. This suggests that a cloud- 
producing scheme needs to be constructed in a framework that also takes the turbulence into consideration. 


1. Introduction 

Cloud is one of the major uncertain quantities in de- 
termining the atmospheric greenhouse warming (e.g., 
Cess et al. 1990; Yao and Del Genio 1999; Bony and 
Dufresne 2005) owing to its role of regulating the radi- 
ative heating through its albedo effect, which tends to 
cool the atmosphere-surface, and through its green- 
house effect, which traps the longwave (LW) radiation 
and tends to warm the atmosphere-surface. The diffi- 
culties in the cloud simulations are not only due to 
coarse model resolutions and imperfect treatments of 
cloud physics but also its complicated interactions with 
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other physics, especially turbulence (e.g., Grenier and 
Bretherton 2001). 

Turbulence effectively transfers heat, moisture, and 
momentum between the surface and the lower atmosphere, 
and modifies the atmospheric stability and vertical distri- 
butions of these quantities. Moist convection and cloud 
formation are strongly influenced by this process. The 
treatments of cloud-topped boundary layers have signif- 
icant effects on the simulations of marine stratocumulus 
(e.g., Park and Bretherton 2009; Lock et al. 2000). 

In earlier GCMs, dry convective adjustment was ap- 
plied to perform the role of turbulence (e.g., Hansen 
et al. 1983). However, it is apparently overly simplified 
and, as shown later, cannot fully substitute for an explicit 
turbulence parameterization. 

In the earlier attempts of parameterizing the turbu- 
lence, the fluxes of heat, moisture, and momentum are 
down the local gradient of the respective mean field. But 
this “local” approach still cannot sufficiently transfer 
these quantities. To more effectively transfer these quan- 
tities vertically, a diffusion term and a countergradient 
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transport term were included in later parameterizations 
of turbulence (e.g., Troen and Mahrt 1986; Holtslag and 
Moeng 1991; Holtslag and Boville 1993). It has been 
emphasized that the nonlocality is included in both the 
proposed diffusivity, via the surface forcing and the 
boundary layer depth, and the additional countergradient 
term (Large et al. 1994, p. 386). 

In this study, we use the Goddard Institute for Space 
Studies (GISS) Model E2 GCM employed in the In- 
tergovernmental Panel on Climate Change (IPCC) Fifth 
Assessment Report (AR5) (G. A. Schmidt et al. 2012, 
unpublished manuscript) to systematically identify the 
influence of turbulence parameterization on the cloud 
simulations. We also present a newer version of the 
turbulence parameterization being developed at GISS 
and compare the cloud and radiation climatologies 
obtained with those obtained from the AR5 version. 

In section 2, we describe the models used. We show 
the experiments conducted in section 3 and the results 
in section 4. A summary and conclusions are given 
in section 5. 

2. The models 

Except for the parameters for computing the thresh- 
old relative humidity RH (t/ 00 ), a minimum RH for 
cloud formation, and the limiting autoconversion rate 
associated with upward motion, we use the same GISS 
Model E2 GCM employed in the AR5 experiments with 
climatological (monthly varying) sea surface tempera- 
tures (SSTs) and sea ice coverage (Rayner et al. 2003) 
averaged from 1975 to 1984. The model has 40 layers in 
the vertical and 2° (latitude) X 2.5° (longitude) hori- 
zontal resolution; the physical time step for cloud and 
turbulence computations is 30 min. 

a. Cloud and moist convection parameterization 

The clouds in the model include stratiform clouds and 
convective clouds. The cloud and moist convection pa- 
rameterization used in AR5 is an improvement on the 
version used in the IPCC Fourth Assessment Report 
(AR4) (Schmidt et al. 2006). 

Numerous improvements have been implemented in 
the AR5 version. For the moist convection (MC), we 
calculate entrainment rates and cumulus updraft speeds 
interactively (Del Genio et al. 2007). MC transports the 
convective condensate vertically. Downdrafts can orig- 
inate from multiple layers and are allowed to detrain at 
lower layers, and drop below the cloud base only if they 
are negatively buoyant. The downdrafts also include 
entrainment and detrainment. The cumulus momentum 
mixing is affected by the cumulus pressure gradient. The 
downdraft mass fluxes are used to calculate a gustiness 


correction to surface fluxes. The convective cloud cover 
depends on the cumulus updraft speed. The rain evapo- 
ration in downdrafts is limited to the extent of saturating 
the downdraft, and one hour is used for the convective 
adjustment time among other minor changes. 

For the stratiform clouds, instead of specifying the 
U QO , we calculate it and use two parameters — U a and 
U b — to achieve radiative balance and a better cloud 
climatology. We set U 00 = U a l[U a + (1 — U a )\p/35], 
where A p is the layer pressure depth and U 00 is fur- 
ther reduced if there is upward motion. However, in 
the moist convective area and in the PBL, we set U OQ = 
1 - 2 U b [<r/q*(T, /?)], where a represents the subgrid 
variability of humidity (Siebesma et al. 2003), q* is the 
saturation specific humidity, and T is temperature; U a 
primarily affects the higher clouds above the 850-mb 
(hPa) level and U b mainly affects the lower clouds below 
the 850-mb level. The ice phase is maintained once 
a cloud forms if T < 0°C and the liquid cloud phase is 
maintained if T > — 4°C. Partial evaporation of cloud 
water in unfavorable conditions for cloud formation is 
allowed, and the optical thickness of precipitation is 
included. Averaged vertical velocity over one physical 
time step is used to modify the autoconversion rate. Ice 
cloud forms if T < -35°C, and the Karcher and Loh- 
mann formula (Karcher and Lohmann 2002) for com- 
puting the saturation vapor pressure below -35°C is 
used among other minor changes. 

b. Turbulence parameterization 

The turbulence parameterization for AR5 experi- 
ments is similar to the version used in AR4 experiments 
(Schmidt et al. 2006). It includes 1) a nonlocal vertical 
transport scheme for virtual potential temperature, 
specific humidity, and other scalars following Holtslag 
and Moeng (1991), inside the PBL and for both stable 
and unstable stratifications; 2) diagnosis of the nonlocal 
vertical profile of the turbulent kinetic energy (TKE) 
according to the large eddy simulation (LES) studies of 
Moeng and Sullivan (1994) that include the PBL height 
in the computation, inside the PBL and again for both 
stable and unstable stratifications; 3) above the PBL, 
employment of the second-order closure (SOC) scheme 
of Cheng et al. (2002) that improves several aspects of 
the standard Mellor-Yamada model (Mellor and Yamada 
1982) and in particular, allows mixing at weak turbulence 
levels; 4) formulation of the turbulence length scale sim- 
ilar to Holtslag and Boville (1993); 5) calculation of the 
PBL height using the TKE criterion (PBL top is where 
TKE decreases to 10% of its surface value), and 6) be- 
tween the surface and the first GCM layer, integrating 
the SOC model equations over the subgrid levels, 
instead of applying the usual interpolation scheme, to 
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find the surface velocity, temperature, humidity, among 
others, which is unique among the GCMs, and is conve- 
nient to add more physics and allows coarse vertical res- 
olutions near the surface. 

Using C as a general variable for the mean virtual 
potential temperature, specific humidity, or winds, and 
c for the respective fluctuating parts, the general equa- 
tion for computing its vertical flux is 


where w 2 is twice the vertical component of the TKE 
e , ( wc) 0 is the value of Wc at the surface, t = 2 els is 
the turbulence time scale with s the dissipation rate 
of e , h is the PBL height, u * is the friction velocity, 
and w* is the Deardorff buoyancy velocity scale that 
also depends on h, defined as 

W* = [ga(we~Xh] m , (4) 


WC = ~ Kc Jz + ( wc ^s’ (1) 

where w is the fluctuating part of the vertical velocity, 
z is the height, and K c is the eddy diffusivity. The first 
term on the rhs of Eq. (1) contains both the local and 
nonlocal contributions, and the second term is the non- 
local countergradient contribution to the vertical turbu- 
lent flux. Note that K c depends on characteristic bulk 
properties of the PBL rather than just local properties 
(see, e.g., Troen and Mahrt 1986; Holtslag and Boville 
1993). 

To further see the influence of the turbulence pa- 
rameterizations on the cloud simulations, we conduct 
experiments using a modified version of the turbulence 
parameterization being developed at GISS. Both the 
turbulence scheme in the AR5 and the newer scheme 
employ the nonlocal model as in Eq. (1) above, but the 
specific forms of the diffusivities and the countergradient 
terms, as well as the related turbulence length scales and 
the PBL height, a critical input to the nonlocal scheme, 
are quite different. This newer scheme differs from the 
AR5 version as follows: 

1) Nonlocal vertical transport scheme for velocity, virtual 
potential temperature, specific humidity, and other 
scalars. Both the AR5 turbulence scheme and the 
newer scheme employ expression (1), but the diffusiv- 
ities K c and the counter gradient flux term (wc) cg are 
modeled differently. 

The AR5 scheme follows Holtslag and Moeng 
(1991). The heat diffusivity, K c , and the counter- 
gradient flux, (Wc) , are derived from the dynamic 
equation for wc, with the third-moment terms 
(related to the transport and the pressure) param- 
eterized using the LES data; the results are 

K c = 0.045vUt, 

and (2) 


where g is the gravitational acceleration, a is the 
thermal expansion coefficient at the surface, and 6 V is 
the turbulent fluctuation of the virtual potential tem- 
perature. 

The newer scheme follows Holtslag and Boville 
(1993), where 


K C 

i^) cg 


KW 


1.2K, 


'K* - 1) 2 ' 

W*(wc) 0 


C w m h 


(5) 


where k = 0.4 is the von Karman constant, w m is 
a velocity scale contributed by both the wind shear 
and the buoyancy, 

w m = (ul + 0.6n4) 1/3 , (6) 


and w t = w m /Pr is the velocity scale for mixing of heat 
or passive scalars, where Pr is the turbulent Prandtl 
number; both w m and w t are characteristic of the 
surface forcing that is exerted to the whole PBL. 

2) Formulation of the turbulence length scale £, defined 
in 


(2e) 312 


(7) 


where Bi = 19.3 is a constant determined in Cheng 
et al. (2002). 

The AR5 scheme follows Holtslag and Boville (1993), 
where 



£ s = kz, and £ 0 = 30 + 270 exp(l - z/1000) (m), 

(9) 

where i s is the length scale near to the surface and is 
the turbulence length scale away from the surface. 
The newer scheme shows that for z ^ h, we use (8) 
while generalizing (9) as 


(3) 


£ 0 = 30 + max(0.3/z - 30,0) exp(l - z/h) (m); 

( 10 ) 


(wc) cg = 0-091 ^-(wc) 0 
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for z<h, we use the formula obtained from the large 
eddy simulation data by Nakanishi (2001), as follows: 


1 

£ 


1 1 1 

p "I" / + / ’ 

i s i 0 i b 


(ii) 


where £ s , £ 0 , and i b are the modified surface length 
scale, the length scale away from the surface, and the 
Deardorff buoyancy length scale, respectively; the 
expressions are formulated to match the LES data, 
and are lengthy so we refer the interested readers 
to Nakanishi (2001) for more details. The buoyancy 
length scale has a dependence on the TKE; the latter 
is parameterized using the LES data (Moeng and 
Sullivan 1994). The TKE is also used in the formu- 
lation of the diffusivities in the local SOC model in 
both the AR5 scheme and in the newer scheme 
(Cheng et al. 2002). 

3) Calculation of the PBL height. 

The AR5 scheme uses the “TKE criterion,” by which 
the PBL top is where the TKE decreases to 10% of 
its surface value. 

The newer scheme uses the “Richardson number 
criterion” according to Troen and Mahrt (1986) 
and Holtslag and Boville (1993), in which the PBL 
height h is determined by 


_ Ri cr [t/(U 2 + V(hf } 

ga[@ v (h) - 0J 


( 12 ) 


where Ri cr (=0.5) is a critical bulk Richardson num- 
ber, U(h) and V(h ) are the horizontal mean wind 
components, 0„ is the mean virtual potential tem- 
perature, and 0 5 is a temperature near the surface that 
may relate to the surface virtual heat flux ( [wd v ) Q de- 
pending on the stability condition, defined as 


0. 


k(z,) + 8.5 


C w6 v)o 




& v( z s)’ 


(wd v ) 0 > o 

K), - °. 


(13) 


where =10 m is the surface layer height. For more 
details, see appendix A of Holtslag and Boville 
(1993). We will further discuss the impacts of the new 
turbulence scheme and the physical reasons behind 
them in section 4. 

In this work we will not include explicit treatments of 
PBL-top entrainment and radiative-cooling-driven tur- 
bulence that are being developed at GISS, but we have 
cloud-top entrainment in the parameterization of large- 
scale clouds (Del Genio et al. 1996). Explicit treatments 
of PBL-top entrainment and radiative cooling effect will 


be included in further development of the turbulence 
parameterization. Despite this, as shown in the results, 
the newer version of the turbulence parameterization 
significantly improves the simulations of cloud and ra- 
diation. 

In recent turbulence parameterizations, the moist 
conserved variables, the liquid water potential temper- 
ature ©/, and the total water mixing ratio Q t are used for 
treating both dry and cloudy layers (Lock et al. 2000; 
Grenier and Bretherton 2001; Bretherton and Park 
2009). But the GISS ModelE uses the potential temper- 
ature and the specific humidity as prognostic variables, 
and the computations of moist convection and large-scale 
clouds are handled separately from the turbulence pa- 
rameterization. Therefore, for consistency and simplicity, 
in this study we use virtual potential temperature and 
specific humidity as prognostic variables, while we are 
working toward using ©/ and Q t in the model. 

3. Experiments 

We use the GISS Model E2 GCM described above 
as the control and then conduct several experiments 
with changes in the parameterization of turbulence, as 
described below: 

1) Experiment CONTROL: the control run as described 
above. 

2) Experiment NO-TURB: an experiment without any 
parameterization of turbulence. 

3) Experiment DRY-CONV: an experiment with dry 
convection to represent the turbulence process. In 
the process of dry convection, vertical mixing of heat, 
moisture, and momentum will occur to restore the 
static stability if the temperature profile is static 
unstable. 

4) Experiment LOCAL: an experiment with only the 
local parameterization of turbulence. 

5) Experiment NEW-TURB: an experiment that in- 
cludes both local and nonlocal parameterizations of 
turbulence, but with certain changes from CONTROL 
as described in section 2. 

We integrate CONTROL and NEW-TURB for 5 yr but 
only 1 yr for NO-TURB, DRY-CONV, and LOCAL, 
since their effects are clearly seen in one year of in- 
tegration and they are not intended for comparison of 
climatologies. CONTROL and NEW-TURB are in ra- 
diative balance by properly choosing the parameters U a 
and U b : U a = 0.7 in CONTROL and 0.72 in NEW-TURB 
and U b = 1.96 in both CONTROL and NEW-TURB. 
NO-TURB, DRY-CONV, and LOCAL use the same 
values for U a and U b as in CONTROL and are not tuned 
to radiative balance. 
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Table 1. Selected global mean and annual mean climate variables, and U a and U b . The planetary albedo and cloud cover are in percent, 
radiative fluxes are in watts per meter squared, and the precipitation is in millimeters per day. 



CONTROL 

NO-TURB 

DRY-CONV 

LOCAL 

NEW-TURB 

Planetary albedo 

29.610 

26.182 

25.680 

27.554 

30.476 

Solar absorption TO A 

240.54 

252.25 

253.97 

247.57 

237.58 

Net LW TOA 

-240.27 

-241.64 

-240.38 

-238.00 

-237.55 

Net heat at surface 

0.27 

10.84 

13.79 

9.79 

0.06 

Precipitation 

3.171 

2.643 

2.925 

3.066 

3.099 

Total cloud 

60.54 

57.63 

55.06 

59.99 

64.96 

High cloud 

31.0 

31.6 

32.5 

34.4 

31.8 

Middle cloud 

15.4 

16.2 

15.8 

16.8 

16.1 

Low cloud 

43.0 

39.9 

32.7 

39.5 

46.9 

U a 

0.70 

0.70 

0.70 

0.70 

0.72 

u b 

1.96 

1.96 

1.96 

1.96 

1.96 


4. Results and analyses 

Here we will compare and discuss the results ob- 
tained from each experiment described in section 3 
with those obtained from CONTROL. Table 1 shows 


some climate variables from these experiments. Note 
that NO-TURB, DRY-CONV, and LOCAL are not 
tuned to radiative balance, and they have very low 
planetary albedo partially because of smaller low cloud 
amounts. 


RH % ANN CONTROL RH % ANN NO-TURB-CONTROL RH % ANN DRY-CONV-CONTROL 



0 10 20 30 40 50 00 70 80 90 100 - 33.9 -7 -5 -3 -1 1 3 5 7 9 14 3 -18 1 -7 -5 -3 -1 1 3 5 7 9 14.1 

Fig. 1. Annual mean and zonal mean of RH (%) of (top) CONTROL, NO-TURB-CONTROL, and DRY-CONV-CONTROL and 
(bottom) NEW-TURB, LOCAL-CONTROL, and NEW-TURB-CONTROL. 
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TOTAL CLOUD % ANN CONTROL 


900 


-90 


ANN NO-TURB-CONTROL 

0| 1 1 l 



ANN DRY-CONV-CONTROL 




TOTAL CLOUD % ANN NEW-TURB 



ANN NEW-TURB-CONTROL 
Oi 1 r 



' 0 % 


-258 


-14 - 6-2 


6 10 14 


Fig. 2. As in Fig. 1, but for total cloud cover (%). 


a. Experiments identifying the turbulent effects 

The experiment without any turbulence parameteri- 
zation has RH and low cloud cover peaking right above 
the surface (Figs. 1 and 2), which is highly unrealistic 
since RH peaks at the top of a well-mixed PBL (Holtslag 
and Boville 1993). The temperature structure in the 
lower atmosphere is statically unstable (not shown), 
suggesting some vertical mixing processes are needed. 

When the dry convection is employed, the moisture 
is sharply reduced near the surface but the strong 
vertical transport is mainly near 900 mb and below 
(Fig. 3). As a result, the humidity and low cloud cover 
below 800 mb are still smaller than those of CONTROL, 
although the vertical distribution of RH and low cloud 
cover are clearly better than the case of NO-TURB 
(Figs. 1 and 2). 

When the local parameterization of turbulence is ap- 
plied, the vertical transports of humidity are still in- 
sufficient above 900 mb in the tropics compared with 
CONTROL, but they are similar in magnitude in the 


middle and higher latitudes (Fig. 3). Consequently, the 
RH and low cloud cover also show similar changes in the 
vertical structures (Figs. 1 and 2). 

The larger RH and cloud cover in the upper tropo- 
sphere in the tropics in NO-TURB and DRY-CONV 
are in part due to weaker deep moist convections that 
dry the atmosphere to a lesser degree because of weaker 
subgrid scale compensating subsidence warming and 
drying in these cases although the shallow convections 
are stronger (not shown), since the atmospheric struc- 
ture is less stable in the lower atmosphere because of less 
turbulent transports. 

b. Comparison of the newer turbulence scheme 
with the control 

We now compare the climatologies of cloud and ra- 
diation obtained from NEW-TURB and CONTROL 
in greater detail. We compare some global mean and an- 
nual mean quantities in Table 1. These values look sim- 
ilar between NEW-TURB and CONTROL, although 
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TURB CH OF LAT HEAT ANN CONTROL ANN NO -TURB- CONTROL ANN DRY-CONV- CONTROL 



TURB CH OF LAT HEAT ANN NEW-TURB ANN LOCAL-CONTROL ANN NEW -TURB -CONTROL 



-100 -60 -20 0 20 40 60 80210 2 -121 -70 -30 10 30 50 70 588.2 - 589 -70 -30 10 30 50 70 299.8 

Fig. 3. As in Fig. 1, but for turbulent changes of latent heat (10 -2 W mb -1 m -2 ). 


NEW-TURB has a higher planetary albedo because of a 
higher low cloud amount. 

Figure 4 compares their geographical distributions of 
total cloud with observations from the International 
Satellite Cloud Climatology Project (ISCCP) (Rossow 
et al. 2002) for January and July. Generally speaking, 
NEW-TURB agrees with observations better, especially 
over subtropical eastern oceans, where fewer clouds 
are simulated in most climate models (e.g., Weare 
et al.1995; Hannay et al. 2009). Further analyses indicate 
the differences of the total cloud are mainly due to low 
cloud simulation (not shown). High clouds associated 
with deep convections may contribute to the differences 
in the tropics (not shown). The larger cloud amounts in 
higher latitudes and smaller cloud amounts in tropical 
oceans obtained by NEW-TURB are closely related to 
the RH in the PBL (Fig. 5). Note that Fig. 5 is for the 
850-mb level and that over the subtropical eastern 
oceans, the RH at 700 mb (not shown) is a better level 
to compare. This is plausibly related to the vertical 


turbulent fluxes of moisture and the diagnosed PBL 
height (see more discussion below). But the shallow 
moist convection (Fig. 6) may also play an important 
role, notably over the subtropical eastern oceans. The 
strength of shallow convection is closely related to the 
stability in the lower PBL, which is affected by the tur- 
bulent fluxes. The increased shallow convective cloud 
cover is small (not shown), but it moistens the atmo- 
sphere at the convective cloud-top level and helps the 
formation of stratiform clouds there. The shallow con- 
vection also intensifies the vertical turbulent transport of 
heat, moisture, and momentum, and as a result deepens 
the cloudy boundary layer (Tiedtke et al. 1988). Note 
that we do not have a separate parameterization for the 
shallow convection as in some GCMs (e.g., Park and 
Bretherton 2009). 

The differences in cloud simulations between NEW- 
TURB and CONTROL clearly will have a large impact 
on the solar absorption (Fig. 7). In general, NEW-TURB 
produces better solar absorption, especially over the 




4970 


JOURNAL OF CLIMATE 


Volume 25 


TOTAL CLOUD % JAN NEW-TURB-ISCCP-0.92 TOTAL CLOUD % JAN CONTROL-ISCCP -4.69 TOTAL CLOUD JAN NEW-TURB-C0NTR0L3.97 



-55 -45 -35 -25 -1ft -5 5 15 25 35 45 66 6 -70.9 - 35 - 25 -15 -5 5 15 25 35 45 64 3 -55 -45 -35 -25 -15 -5 5 15 25 35 45 55 


Fig. 4. Total cloud cover (%) of NEW-TURB-ISCCP observation, CONTROL-ISCCP observation, and NEW-TURB -CONTROL for 

(top) January and (bottom) July. 


subtropical eastern oceans and the southern oceans. 
But the solar absorption over the subtropical eastern 
oceans is still too large compared with the observation. 
More complete treatments of cloud-top-driven turbu- 
lence may be needed to reduce this quantity. The solar 
absorption is generally too strong over the southern 
oceans in climate simulations, which may affect the 
cloud feedback in GCMs (Trenberth and Fasullo 
2010). The solar absorption over the Arctic oceans in 
July is apparently too large in both NEW-TURB and 
CONTROL. This excessive solar absorption helps 
melt the sea ice in July and is designed to prevent ex- 
cessive growth of sea ice in the coupled runs in which 
the sea surface temperatures are predicted. However, 


this is only a temporary fix until a better sea ice 
transport scheme and stronger poleward heat trans- 
ports can avoid the excessive sea ice growth in coupled 
runs. 

To find out the differences in the turbulent processes 
between NEW_TURB and CONTROL that can, at 
least qualitatively, explain the differences in the simu- 
lations of cloud and radiation, we show the boundary 
layer heights h in Fig. 8, in which h obtained in NEW- 
TURB is higher over land and over oceans in middle and 
higher latitudes but is lower over oceans in the tropics. 
However, h is higher over subtropical eastern oceans in 
NEW-TURB than CONTROL. So, except for over 
land, the differences in h correlate with the differences 


RH 850MB % JAN NEW-TURB 64 83 



10 20 30 40 50 60 70 00 90 101.5 


RH 850MB % JAN CONTROL 


67.04 RH 850MB JAN NEW-TURB -CONTROL -2.21 



Fig. 5. 850-mb RH (%) of NEW-TURB, CONTROL, and NEW-TURB-CONTROL for (top) January and (bottom) July. 
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— 1— I 1 1 II 


-30 -26 -16 -6 


Fig. 6. As in Fig. 5, but for shallow convective frequency (%). 


SHLW CONV FREQ % JAN CONTROL 41.73 SHLW CONY FREQ NEW-TURB -CONTROL 8. 31 


SHLW CONV FREQ % JUL CONTROL 42.05 SHLW CONV FREQ NEW -TURB- CONTROL 7 95 


SHLW CONV FREQ % JAN NEW-TURB 50.04 


SHLW CONV FREQ % JUL NEW-TURB 50.00 


in cloud simulations quite well. Higher h tends to 
transfer moisture higher and produces higher RH and 
cloud cover. The differences in h are clearly affected by 
the different methods of determining the PBL top be- 
tween NEW-TURB and CONTROL. A test experiment 
with NEW-TURB using the method of calculating the 
PBL top as in CONTROL further confirms this point 
(not shown). Nevertheless, larger h tends to produces 
larger turbulent fluxes [see Eqs. (l)-(5)] and has sig- 
nificant effects on RH and the low cloud amount. Over 
land, this relationship is not clear because of compli- 
cated interactions with the surface physics. The bound- 
ary heights over midlatitude and high-latitude oceans 
seem to be on the high side, but they appear to be similar 


in magnitude to the observations of von Engeln and 
Teixeira (2012, manuscript submitted to J. Climate). The 
other two main differences in nonlocal vertical transport 
and length-scale formulation between NEW-TURB and 
CONTROL appear to have only produced small im- 
provements in the results (not shown). 

Since the precipitation is associated with the cloud 
processes, we compare the precipitation climatologies 
with observations of the Global Precipitation Clima- 
tology Project (GPCP) (Huffman et al. 1997) in Fig. 9. 
Both NEW-TURB and CONTROL produce the pre- 
cipitation reasonably well, although NEW-TURB pro- 
duces weaker precipitation than what is observed over 
the Amazon area in July. 
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Fig. 7. Solar absorption at top of atmosphere (TOA) (W m -2 ) of NEW-TURB minus the Clouds and the Earth’s Radiant Energy System 
(CERES) observation, CONTROL-CERES observation, and NEW-TURB-CONTROL for (top) January and (bottom) July. 
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Fig. 8. As in Fig. 5, but for PBL height (m). 


The latitude-pressure cross sections of differences in 
RH, total cloud cover, and turbulent latent heat be- 
tween NEW-TURB and CONTROL are shown in Figs. 
1-3. NEW-TURB vertically transports heat and mois- 
ture not as deep as in CONTROL in the tropics, but 
deeper in mid- and higher latitudes. This is likely related 
to the PBL height calculated. The RH and low cloud 
cover have vertical distributions reflecting these differ- 
ences in the turbulent transports. This may also con- 
tribute to the differences in their cloud and radiation 
climatologies. 

We also note that NEW-TURB produces more re- 
alistic diurnal variations of convective precipitation. 
Over land in the tropics, the convective precipitation 


peaks around 1500 LT instead of around 1200 LT in the 
AR5 version (not shown). Better interactions between 
the moist convection and PBL are important for simu- 
lating the diurnal variations of convective precipitation 
(Rio et al. 2009). Further analyses are needed for this 
aspect of cloud simulation and will be a subject for fu- 
ture studies. 

5. Summary and conclusions 

Using the GISS Model E2 GCM employed in AR5 
experiments, we have carried out a series of experiments 
to explore the impact of turbulence parameteriza- 
tions on the cloud simulations. Without the turbulence 
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Fig. 9. Precipitation (mm day x ) of NEW-TURB, GPCP observation, and NEW-TURB-CONTROL for (top) January 

and (bottom) July. 
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parameterization, the RH and the low cloud cover peak 
unrealistically close to the surface; with the dry con- 
vection parameterization, these two quantities improve 
their vertical structures, but the vertical transport of 
water vapor is still weak in the PBL; with only local 
turbulence parameterization, the effects are similar to 
the case of the dry convection in the tropics; with both 
local and nonlocal turbulence parameterizations, the 
RH and low cloud cover have better vertical structures 
in the tropics due to more significant transport of water 
vapor in the PBL. This suggests that a cloud-producing 
scheme needs to be constructed in a framework that also 
takes turbulence into consideration. 

The climatologies of cloud and radiation obtained 
from a newer version of turbulence parameterization 
are also compared with those obtained from the AR5 
version. This newer version of parameterization has 
significant improvements in cloud and radiation sim- 
ulations, especially over the subtropical eastern oceans 
and the southern oceans. The differences of the di- 
agnosed PBL height between the newer version and 
the AR5 version appear to correlate well with the 
differences of the total cloud distribution over oceans. 

Although this newer version of turbulence parame- 
terization produces better simulations of cloud and ra- 
diation, its vertical transports of heat and moisture in the 
tropics are not as deep as in the AR5 version. 
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